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Abstract The Monge-Kantorovich mass transfer problem is equivalently formulated as an 
optimal control prblem for the mass transport equation. The equivalency of the two problems 
is establish using the Lax-Hopf formula and the optimal control theory arguments. Also, it 
is shown that the optimal solution to the equivalent control problem is given in a gradient 
form in terms of the potential solution to the Monge-Kantorovich problem. It turns out that 
the control formulation is a dual formulation of the Kantrovich distance problem via the 
Hamilton- Jacobi equations. 



1 Introduction 

Monge mass transfer problem is that given two probability density functions po(x) > and 
Pi{%) > of x G R d , find a coordinate map M such that 



(1.1) / pAx) dx = / p (x)dx 

J A JM(x)eA 

for all bounded subset A in R n . If M is a smooth one-to-one map, then it is equivalent to 

(1.2) det(VM)(x)pi(M(x)) = p (x) 

where det denotes the determinant of Jacobian matrix of the map M. Clearly, this problem is 
underdetermined and it is natural to formulate a costfuctional for the optimal mass transfer. 
The so-called Kantorovich (or Wasserstain) distance between p and p\ is defined by 



1.3) dfapj-mtf c(x-M(x)) Po (x)dx 

JR d 



where c is a convex function and c(x — y) = c(\x — y\) with c(0) = 0. For example c(x — y) = 
~\x— y\ p is for the LP Monge-Kantorovich problem (MKP). Whenever the infimum is attained 
by some map M, we say that M is an optimal transfer for the Monge-Kantorovich problem. 
The Kantorovich distance is the least action that is necessary to transfer p into p\. 

The mass transport problems have attracted a lot of attentions in recent years and have 
found applications in many fields of mathematics such as statistics and fluid mechanics 
(e.g., see [21 [5] and [11] for extensive references). From a more scientific point of view the 
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Kantorovich distance provides a valuable quantitative informations to compare two different 
density functions and it has been used in various fields of applications [3]. 
It is shown e.g., in [H HI [HI [5] that the optimal map M is given by 

(1.4) Dc(x - M(x)) = Vm(x) 

for a potential function u, where Dc denotes the derivative of c. In fact u is the optimal 
solution to the Kantorovich dual problem (2.2). If c is uniformly convex, then we can solve 

(1.4) for M in terms of Vm. For example for LP MKP 

x - Mix) = \Vu(x)\ q ~ 2 Vu(x) with - + - = 1. 

P Q 

For L 2 MKP, it follows from (1.2) and (1.4) that if ip = &f-u(x), then M(x) = x-Vu = Vip 
and thus ifi satisfies the Monge- Ampere equation 

(1.5) det(HiP)(x) Pl (Vip) =po(x), 

where Hip is the Hessian of ip. 

In [3j the L 2 MKP is equivalently reformulated as an optimal control problem: 



-/ [ p(t,x) \V(t,x)\ 2 dxdt over vector field V = V(t, x) 

2 JO JR d 



(1.6) d(p ,pi) = min 
subject to 

Pt + V- (pV) =0 

(1.7) 

p(0, x) = P o(x) and p(l, x) — pi(x) 

Moreover if V(t, x) is an optimal solution to (1.6)— (1.7) and the Lagrange coordinate X(t; x) 
satisfies 

^-X = V(t,X(t;x)), X(0;x)=x, 
at 

then M(x) = X(T; x). 

The contribution of this paper is that we will show that the optimal vector field V to 
problem (1.6)-(1.7) is given by 

(1.8) V(x,t) = V x $(t,x) 

where the potential function satisfies the Hamilton- Jacobi equation 

(1.9) ^ + I|V0| 2 = O, $(Q,x) = -u(x) 

and u determines the optimal map M in (1.4). Thus, (1.8)-(1.9) is an optimal feedback 
solution to control problem (1.6)— (1.7), i.e., given p , p\ first we determine ip by (1.5) and 

I oc I ^ 

let u = ip and then determine V by (1.8)— (1.9). 
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Moreover, it will be shown that 



d{po,P\) = min 

subject to 



(pi{x)v(x) — po(x)(j)(0, x)) dx over v 

R d 



0t + ^|V0| 2 = O, <j)(l,x) = v(x). 

It is the other control formulation of the L 2 MKP and is an optimization problem over the 
potential fuction v subject to the Hamilton- Jacobi equation. 

For the non-quadratic c case, the (generalized) optimal control problem is formulated as 

(1.10) min / / pit, x) c(V(t, x)) dxdt 

Jo JR d 

subject to (1.7). In this case the optimal vector field V is given by 

(1.11) V = De*(V(j>) 
where c* is the convex conjugate function of c defined by 

c*(y) = sup {x ■ y — c(x)}. 



For U> MKP 



and 



c(x) = -\x\p, xeR d , c*(y) = -\y\«, yeR d , 
V Q 



V(t,x) = \V x 0(t,x)r 2 V x (f>(t, 



x) 



where — I — = 1 and pG (1, oo). The potential function = <f>(t, x) satisfies 
V Q 

(1.12) 4 + c*(V0)=O, 0(0, ac) = -u(x). 
If c is uniformly convex, then from (1.4) 

M{x) = x- Dc*(Vu(x)). 

Thus, from (1.2) u satisfies 

(1.13) det(VM)(i)pi(a; - Dc*(Vu{x))) = p (x). 

For L 2 MKP (1.13) is reduced to (1.5). Hence the optimal solution to (1.10) subject to (1.7) 
is given in the feedback form (1.11)-(1.13). 

An outline of our presentation is as follows. In Section 2 the basic theoretical results 
concerning the MKP problem is reviewed following [5]. Then equivalent variational formu- 
lations (2.5) and (2.8) for the potential function are then derived using the duality and the 
Lax-Hoph formula. In Section 3 we present formal arguments that show the feedback solu- 
tion (1.7)— (1.9) to (1.5)— (1.6). In Section 4 we validate the steps in Section 3 mathmatically 
for L 2 MKP. In Section 5 we present the proofs for the general case. 
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2 Variational Formulations 



In order to present our treatment of the MKP problem, we first recall a basic theoretical 
result in this section. The following relaxed problem of (1.3) is introduced by Kantorovich. 
Let M. be a class of random probability measures p on R d x R d satisfying proj^ = p dx 
and projz/i = p\ dy. Then we define the relaxed cost-functional 



(2.1) J(fi) = / c(x — y) dfi(x,y) over M.. 

JR d xR d 

Consider the dual problem of (2.1); maximize 
(2.2) 



u(x)po(x) dx + / v(x)p\(y)dy 

R d JR d 



subject to u{x) + v (y) < c(x — y). 

The point of course is that the Lagrange multiplier associated with the inequality in (2.2) 
solves problem (2.1). The following theorem pfl HI [8] provides the solution to (2.2) and 
(1.3). 



Theorem 2.1 

(1) there exists a maximizer (u,v) of problem (2.2). 

(2) (u,v) are dual c-conjugate functions, i.e., 

u{x) = inf (c(x — y) — v{y)) 

v(y) = inf (c(x — y) — u(x)) 

(3) M(x) satisfying Dc(x — M(x)) = Vu(i) solves MKP problem. 
It follows from Theorem 2.1 that (2.3) is reduced to maximizing 

(2.4) J(u)=l u(x)p (x)dx+ v(y)pi(y)dy 

JR d JR d 

over functions u, where v is the c-conjugate function of u. The c-conjugate function of a 
function u is defined by 

v (y) = inf (c(x — y) — u(x)). 

It is easy to show that the bi c-conjugate function u of u satisfies u > u a.e. and thus the 
maximizing pair (u,v) of (2.4) is automatically c-conjugate each other. Similarly, we have 
the equivalent problem of maximizing 

(2.5) J(v)= u(x)p (x)dx+ v(y)p 1 (y)dy 

JR d JR d 

4 



where 

u(x) = inf (c(x — y) — v(y)). 
Let c* be the convex conjugate of c, i.e., 

c*(x) = sup((x,y) - c(j/)). 
y 

By the Lax-Hopf formula [6], if <fi is the viscosity solution to 

(2.6) t + c*(V0) = O, <p(l,y) = v(y) 
then 

(2.7) 0(0, x) = sup (f (y) — c(x — y)) = —u(x). 

y 

Thus, Problem (2.2) can be equivalently formulated as maximizing 



(2.8) J(v)= (pi(x)v(x)-po(x)<P(p,x))dx 

subject to (2.6). 

3 Derivation of Optimal Feedback Solution 

The optimality condition of (2.8) subject (2.6) is formally derived as follows. We define the 
Lagrangian 

(3.1) LU,X) = J - f 1 f U t + c*(V<f)))\dxdt. 

JO JR d 

By applying the Lagrange multiplier theory the necessary optimality is given by 



(3.2) 



L^(0, X)(h) = f 1 [ (A* + {Dc\V(j)) \) x )hdxdt 

JO JR d 

- I (h(l,x)\(l,x)—h(0 1 x)\(0,x))dx+ [ (h(l, x)pi(x) - h(0, x)p (x)) dx = 
Jn d Jn d 



for all h G Cq([0, 1] x R d ). Hence the necessary optimality reduces to 

\ t + (Dc*(V4>) X) x = 

(3.3) 

A(0, x) = p (x), A(l, x) = pi(x). 
This implies that if we let V(t,x) = Dc*(V<j)(t,x)) in 

pt + (Vp) x = 0, p(0,x)=p , 



then p(l,x) = pi(x). Moreover, we can argue that 

(3.4) / (p u(x) + pi(x)v(x))dx— / / p(t, x)c(V(t, x)) dxdt = 

JR d Jo J R d 

since 

c(v) = (yJ,Dc*(yJ)) - c*(v x 0). 

It follows from (3.3)-(3.4) that V = Dc*(V(p) is the optimal solution to (1.10) subject 
to (1.6). In fact, for sufficiently smooth pair (p, V) satisfying (1.6), we define Lagrange 
coordinate X(t,x) by 

^-X = V(t,X(t,x)), X(0,x) = x, 

Then for all test function / 

/ / f(t,x)p(t,x)dxdt— / / f(t,X(t,x))p (x)dxdt. 
Jo Jn d Jo JR d 

(3-5) 

f(t,x)p(t,x)V(t,x)dxdt= [ [ V(t,x(t))f(t,X(t,x))p (x)dxdt. 



10 JR d Jo JR d 

Note that (1.6) and (3.5) imply that M(x) = X(l,x) satisfies condition (1.1). Letting 
/ = c(V) in (3.5), we have 

1 [ p(t,x)c(V(t,x))dxdt= f 1 f c(V(t,X(t,x))p (x)dxdt. 

JR d Jo JR d 

f 1 [ c(^-X(t,x))po(x)dxdt> [ c(X(l,x)-X(0,x))po(x)dx. 
o JR d dt J 

where we used the Jessen's inequality. Since M(x) is the optimal solution to (1.3), it follows 
that 

(3.6) I> [ c(x-M(x))p (x)dx = d(p , Pl ). 

jR d 

From (3.4), (3.6) and Theorem 2.1 

(3.7) d(p , pi) = / (p u(x) + pi(x)v(x)) dx = / / p(t, x)c(V(t, x)) dxdt < I. 

jR d Jo JR d 

for all pair (p,V) satisfying (1.6). That is, (p,V) is optimal. 



4 Proof of (3.3)-(3.4) 

In this section we give a proof for the steps of deriving the optimality condition (3.3) and 
equality (3.4) in the case when p — 2, i.e., c(\x — y\) — \ \x — y\ 2 . Suppose v G iy 1,00 (i? d ) 
and v is semi-convex. Then it follows from the Lax-Hopf formula 

<j>{t, x) = sup {v(y) - (1 - t)c{^-)} 
y l — t 
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(e.g., see [6]) that (2.6) has a unique solution <fi G W 1,oo ([0, 1] x R d ) with 
(4.1) \<p{t)\ w i,c* < \v\ w i,oo and \4> t {t)\ L ™ < 



2 



and 

(4.2) cj>(t, x + z) - 2(j>(t, x) + <j){t, x- z)>-C \z\ 2 for all t G [0, 1] and x, z G R d . 

C 

where we assumed v H \x\ 2 is convex. Let r be the solution to (2.6) with r (l) = v + r h 

C 

for h G Co(-R d ). Assume <^ T (1) + "^"l x | 2 be convex for \t\ < 1 and thus (4.2) holds for cjf . 

Step 1 Since y — > c(x — y) — is coersive, for each x E R d there exist y, y T G _R d such 
that 

<f>(0,x) = v(y) - c(x - y), (j) T (0,x) = (v + rh)(y T ) - c(x - y T ). 

Thus 

0(0, x) > v{y T ) — c(x — y T ) = (x + r h){y T ) — c(x — y T ) + r h{y T ) 

and 

0(O,x)-0 r (O,a;) >r% r ) 



Similarly 



r (O,x)-0(O,x) >r%). 



Hence 

(4.3) |0 T (O,-)-0(O,-)|oo<r|/i| oo . 
Step 2 Note that 

(4.4) (0 T - d>) t + i(V0 T + Vcf>) ■ (V0 r - V0) = 0. 
Let rj e , e > be the standard molifier. Then 

(4.5) |V(77 e *^)U < IV0U 
and 

(4.6) V(?y e * 0) -> V</> a.e. as e ^ 0+. 



Moreover (4.2) implies 
Thus 



D 2 (^ *<j))> -C. 
(V(r/ e *0), W) < / ifjdx 

JR d 



for ip G W l,1 (R d ) and ip > a.e. in R d . Thus from (4.5)-(4.6) and the Lebesgue dominated 
convergence theorem, letting e —>■ + 

(4.7) ( V0, V^) < rfC / ^dac 

for -0 G and ^ > a.e. in It now follows from (4.4) and (4.7) that 

j t \4> T {t, ■) - .)|i > -dC |0 T (t, ■) - 0(t, Oli, r (l) = 0(1) + r 

and thus 

(4.8) |^(O,0 -0(0, Oil <re dC \h\ x . 
Since from (2.6) 



we have 



f 1 [ -\V(f) T \ 2 dxdt= [ U T (0,x) -v(x))dx, 

JO JR d 2 jRd 

f 1 [ \V<p T \ 2 dxdt^ f 1 [ \\/(p\ 2 dxdt. 

JO JR d Jo JR d 

as t — > 0. Since L 2 ((0, 1) x R d ) is a Hilbert space, this implies that 
(4.9) f 1 [ | V0 r - V0| 2 dxdf -> 

Jo JR d 

as r — > 0. 

Step 3 For e > let us consider 

A i + (V0A) a; = eAA, A(0) = Po 
Since is Lipschitz on [0, 1] x R d , 

*-> / [e(VA, W) - (V0(V)A,V^)]dx 

defines an integrable, bounded, coersive form on H l (R d ) x H l (R d ) and thus it follows from 
the parabolic equation theory (e.g., see [T21 HO]) that there exits a unique solution X e G 
H\0,l;L 2 (R d ))nL 2 (0,l;H 2 (R d )) provided that p G H 1 (R d ) D L 1 (R d ) D L°° (R d ) . Moreover 

|A e (t)|i < |po|i, 
(4-10) | (!A e (l)|l - L 

|A e (t)|oo < e ct \po 











ID- 






JO 
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For the last estimate we have from (4.7) 

- 4- \k\ p < p — (V0, v|A e n < {p ~ l)dC 

p at p p 

for p > 1 and thus 

I \ I ^dCt I I 

Thus X e is uniformly bounded in L 2 ((0, 1) x R d ). Hence there exists a A G L°°(0, 1; L l (R d ) n 
L°°(R d )) and subsequaence of A e (denoted by the same) such that A e converges weakly to A 
in L 2 ((0, 1) x i? d ) and A e (l) -> A(l) in L 2 (Q). Since for ^ e Cg([0, 1] x 

(4.11) Z 1 / (A e ^ t + A e V0- V^-eVA e - V^)dxdt = [ (p o ip(0) - A e (l)^(l)) dx 

JO JR d JR d 

it follows from (4.10)-(4.11) that letting e -> 0+ 

(4.12) f 1 / (A^ + AV0 • W>) = / (poV'(O) - A(l)^(l)) dr. 

JO JR d JiJ d 

Hence A is a weak solution to 

(4.13) A t + (AV0) = O, A(0) = Po . 

Next we show that (4.13) has the weak unique solution in L°°(0, 1; L l (R d ) nL 1 (R d )). Let 
r] t , e > be the standard molifier and consider the adjoint equation 

(4.14) ^ + Vfo * 0) • Vij = f e C o °°([0, 1] x i? d ), = 0. 
Then, (4.14) has a smooth unique sulution -0 and J = |Vt/>| satisfies 

(4.15) J t + V( V e*(p)- VJ + D 2 ( Ve *4>)J = V/, J(1) = 0. 

Since -0 has compact support, J has a positive maximum over [0, 1] x R d at some point 
(t , x Q ). If < t < 1, then from (4.15) 

Jt(to,x ) < and VJ(t ,^o) = 0. 

Thus, 

L> 2 (7/ e *0)J(i o ,zo) > Vf(to,x ) 
Since from (4.2) -D 2 (?7 e * 0) < — C, this implies 

(4.16) |V^|oo = ^o^o)<^^ 

Let A, A is two weak solutions to (4.13). Then, it folows from (4.12) and (4.14) that 
t I (\-~X)fdxdt= f 1 I (A - A)(V(?7 e * 0) - VcfyVip dxdt. 

JO JR d Jo JR d 

9 



By letting e — > + , it follows from (4.5)-(4.6), (4.16) and the Lebesgue dominated conver- 
gence theorem that 

J J (X-X)fdxdt = Q 

for all / G C °°([0, 1] x R d ) and therfore A = A. 

Now, let A T be the solution to (4.13) associated with (jf . Since A T is uniformly bounded 
in L°°(0, 1; L l {R d ) n L°°(R d )), there exists a A* G L°°(0, 1; L\R d ) n L°°(R d )) such that A r 
converges weakly to A* in L 2 ((0, 1) x R d ) and A r (l) converges weakly to A*(l) in L 2 (R d ) as 
r -> 0. Note that 

(4.17) Z" 1 / (A r ^ + (A r V0 + A r (V0 r -V ( /))-V^)M= / (po^(0) - A r (l)^(l)) dx. 

JO Ji? d JR d 

Since A r is uniformly bounded in L°°(0, 1; L 1 (R d ) (1 L°° (R d )) , it follows from (4.9) and (4.17) 
that A* is the weak solution to (4.13). Since (4.13) has the unique weak solution, we conclude 
A r converges weakly to A in L 2 ((0, 1) x R d ) and A r (l) weakly to A(l) in L 2 (R d ) as r — > 0. 

Step 4 Note that 

(4.18) f 1 [ A r (^ + V0 r • W>) dxdt = [ (^(0)po-^(l)A r (l))dx. 

JO JR d JR d 

for all V e VP 1,oo ((0, 1) x # d ). Since 

_ 0)i + V 0T . y(0 r _ 0) _ 1 | V ^r _ V0 |2 = q 

by setting *0 = r — in (4.18), we obtain 

((0 r (O) - <£(0)Wa;) - t X T (l)h(x)) dx = \ f 1 f X T | V0 r - V0| 2 dxdt 

R d 2 JO JR d 



Similarly, since 



( r _ + V(/) . V(0 . _ 0) + 1 | V0 r _ V 0| 2 = 0, 



we have 

/ ((</> r (0) - 0(O)Wx) - r A(lWx)) dx = -\ f 1 f A | V0 r - V0| 2 dxdt. 

JRd 2 JO JR d 

From (4.3) and (4.8) there exists a subsequence of - — that converges weakly in L 2 (R d ) 

as t — > 0. Since A, A r > a.e. in (0, 1) x i? d and A T (1) converges weakly to A(l) as r — > 0, 
we have 

l im / r (O) -0(0) ,x dx = /" A (i)/ l ( x ) dx . 

Hence 

(4.19) J'(v)(h)= [ (p 1 -X(l))h(x)dx. 

JR d 
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Step 5 Assume v attains the minimum of J(v) in (2.8) and v is Lipshitz and semi-convex. 
Then J'(v)(h) = for all h G Cl(R d ) and thus from (4.19) A(l) = p x a.e., where A is the 
weak solution to (4.13) with = 0. Thus, (3.3) holds with p = A. Since 

t + V0-V0-^|V0| 2 = O, 

it follows from (4.12) with -0 = that 

/ (po(x)u(x) + pi(x)v(x)) dx = [ [ - X(t, x)\V(t, x)\ 2 dxdt. 
jR d Jo JR d 2 

which shows (3.4). 



5 General Case 

In this section we prove (3.3)-(3.4) for the general case c(x — y) = —\x — y\ p . 1 < p < oo. 

P 

Assume v is Lipschitz and semi-convex. It the follows from [5] that 

(5.1) t + I|V0|' ? = O, <P(l,x)=v(x) 

has a unique viscosity solution <fi e W 1,oo ((0, 1) G R d ) satisfying (4.2). 

Step 1 For h G Cq(R ) let v T = v + t h. If for x G R d , y T = y T {x) G R d attains the maximum 
of y — >• v T (y) — c(x — y), then 

x-y T = -\Vv^y T )r 2 Vv T (y T ) 

Since v T G W 1,oc (R d ) , |y r (x) — x\ < a for some a uniformly in x and r. Since as shown in 
Section 4 

r%)<7r T (O,x)-0(O,aO<-T% T ), 



|0 r (O,a;)-0(O,x)|(ix = r / (\h(y(x))\ + \h(y T (x))\) dx. 

R d JR d 

Since h is compactly supported, it follws that there exists a constant M (depends on h) such 
that 



(5.2) / |0 T (O,x) -0(0, x) | dx< Mr. 

JR d 

Since from (5.1) 

f 1 I -\V(j) T \ q dxdt= [ (0 r (O,x) -cf) T (l,x))dx, 
Jo JR d a JR d 



we have 



f I (|V0 T | 9 - \V(j)\ q )dx^0 

Jo JR d 
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as t — > 0. Hence 



t I \\V(t> T \^V(j) T \ 2 dxdt^ f 1 f \\V(j)\^V<p\ 2 dxdt 

JO JR d JO JR d 

asr^O and therefore 

|V0 T | 2 ^V0 r -> |V0|^V0 in L 2 ((0, 1) x 

Moreover, there exists a subsequence (denoted by the same) of r such that V(f) T (x) — > V0(x) 
a.e. in Since by the Lebesgue dominated convergence theorem 

(5.3) \V(f> T \ q - 2 V<p T -> |V0| 9 - 2 V0 in L 2 ((0, 1) x i? d ). 

Step 2 We assume that for p > 2 (i.e., g < 2) | Vv(-) | 2 > c > a.e. in R d . Then it follows 
from (5.1) that |V0(t, -)\ 2 > c a.e. in for t G [0, 1]. Let 

j e = \v( Ve *<p)\ q - 2 v( Ve *<p) 

Then, 

V • J e = | V(7/ e * 0)r 4 (| Vfo * 0)| 2 A(r/ e * 4>) + (q - 2) Vfo * ^'[D 2 ^ * 0)]V(r? e * 0)) 

Since from (4.2) D 2 (j] e * <fi) > — C, there exists a positive constant C q such that V • J t > —C q 
and thus 

(J e ,Vip) < C q [ i)dx 

jR d 

for ip G W 1,1 {R d ) and -0 > a.e. in R d . It thus follows from (4.5)-(4.6) and the Lebesgue 
dominated convergence theorem that 

(5.4) (|V0|^ 2 V0,V^) < C q I i)dx. 

jR d 

for ip G W 1 ' 1 (R d ) and ip > a.e. in R d , by letting e — > + . 
Step 3 Using the same arguments as in Step 3 in Section 4, 

A t + (|V0 r r 2 V0 T A) x = 

has the unique weak solution A r G L°°(0, 1; L\R d ) n L°°(R d )), i.e., 

(5.5) f 1 / A r (^ + |V0 r | 9 " 2 V0 r • V4>)dxdt = [ (i)(<d)p Q -iP(l)\ T (l))dx. 

Jo JR d JR d 

for all -0 G W 1 ' oo ((0, 1) x i? d ). Moreover A r converges weakly to A in L 2 ((0, 1) x R d ) and 
A r (l) converges weakly to A(l) in L 2 (R d ) as r — > 0. Since 

(0 T - 0) t + |V0 r |^ 2 V ( /) r • V(0 r - 0) + h = 

12 



where 

/ x = i|V0 T | 9 - ^|V0|" - \V<j) T \ q - 2 V<j) T • (V0 r - V0) < 
By setting -0 = r — in (4.5), we obtain 

/ ((0 r (O) -0(O))p o (z) -rA r (l)/i(x))da; = - / A r h dx 

JR d JR d 

Similarly, since 

(<p T - 4>)t + |V0r 2 V0 ■ v(0 r - <t>) + h = o 

where 

j 2 = -^|V0 T | 9 - i|V0| 9 - |V0| 9_2 V0 • (V0 T - V0) > 0, 

we have 

J Rd ((<P T (0) ~ m)po(x) - r q{l)h(x)) dx = -J\I 2 dx. 

Since A, A r > a.e. in (0, 1) x R d and A r (l) converges weakly to A(l) as r — > 0, it follws 
that 

lim f - - - — ^^—p (x) dx = [ \(l)h(x) dx. 

Thus 

(5.6) J'( v ){h)= ( ( Pl -X(l))h(x)dx. 

JR d 

Assume v attains the minimum of J(v) in (2.8) and v is Lipshitz and semi-convex. Then, 
J\v){h) = for all h G Cl(R d ) and therrrfore from (5.6) A(l) = p\ where A is the weak 
solution to 

(5.7) A t +(A|V0r 2 V0)* = O, \(0) = Po 
Thus, (3.3) holds with p — A. Since 

t + |V0r 2 V0-V0- V</>l 9 = o, 

V 



it follows from (5.5) with -0 = that 



/ (po(x)u(x) + pAxWx)) dx = f f - \(t,x)\V(t,x)\ p dxdt. 

JR d JO JR d V 



IR d JO JR d p 

which shows (3.4). 
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